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INTRODUCTION 


This report provides simulation results for constant envelope pulse shaped 8 Level 
Phase Shift Keying (8 PSK) modulation for end to end system performance. In order to 
increase bandwidth utilization, pulse shaping is applied to signals before they are 
modulated. This report provides simulation results of power spectra and measurement of 
bit errors produced by pulse shaping in a non-linear channel with Additive White Gaussian 
Noise (AWGN). The pulse shaping filters can placed before(Type B) or after(Type A) 
signals are modulated. Three kinds of baseband filers, 5th order Butterworth, 3rd order 
Bessel and Square-Root Raised Cosine with different BTs or roll off factors, are utilized in 
the simulations. 

The simulations were performed on a Signal Processing Worksystem(SPW). This 
project is conducted at New Mexico State University (NMSU) in the Space 
Communication and Telemetry Laboratory in the Electrical and Computer Engineering 
Department. 


BACKGROUND 


The most bandwidth-efficient communication methods are imperative in coping 
with the congested frequency bands. Pulse shaping methods have excellent effects for 
narrowing the bandwidth and increasing band utilization. For 8 PSK, a non constant 
envelope spectrum shaping method was previously simulated [5]. In these simulations, 
the spectrum of filtered non-constant envelope 8PSK spreads as the signals go through the 
Solid State Power Amplifiers SPA). The position of the baseband filters for the pulse 
shaping are crucial, since if the filter is after the modulator(this will be referred to as type 
A) and before the modulator(this will be referred to as type B) have different effects for 
narrowing the bandwidth and producing bit errors. The constant envelope 8 PSK and pulse 
shaping filters(Type B) are used to perform the simulations on spectrum and bit error rates 
in this project. 
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PROCEDURE 


8 PSK Theoretic Analysis And Pulse Shaping Methods 

This project provides an end to end simulation system. The random binary signal is 
combined via a Gray code. Then the 8-level signal passes through a pulse shaping filter, is 
modulated, passed through a channel and then a receiver. Figure 1 shows the system 
architecture. 



data 


Figure 1 8PSK End To End System Block Diagram 


8 PSK Modulation And Pulse Shaping Filters 

In 8 level Phase-Shift Keying, information is contained in the phase, and one 
symbol has 3 bits. The phase of the carrier can take one of 8 discreet equally spaced 
values, such as 22.5°, 67.5°,1 12.5°, 157.5°, 202.5°, 247.5°, 292.5°, 337.5°. An 8 PSK signal, 
s(t), can be expressed as follows: 


s(t)=A c exp[}(w c t+ 6 )]=A c cos[w c t+ 0 (t)]+j A c sin[w c t+ 0 (t)] 


Where the envelope A c is constant in this project, w c =2 n f c , where f c is the carrier 
frequency in hertz, and 6 is the phase which takes one of the 8 discrete values in degree 
corresponding to the 3 bits Gray code. Table 1 shows the relationship between Gray codes 
and 0 . 
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Binary Gray Codes 

<9 (Degree) 

000 

22.5° 

001 

67.5° 

011 

112.5° 

010 

157.5° 

110 

147.5° 

111 

247.5° 

101 

292.5° 

100 

337.5° 


Table 1 Gray Codes And Phase 


In SPW, the Gray codes and phase 6 can be produced by the design shown in Figure 2. 



Figure 2 Producing Phase Signals Block Diagram 

Where Cl=l, C2=2, and C3=4. The Number to Symbol block which has the function as 
Table 1 completes the conversion from Gray codes to phase in degrees. 

NRZ-L Ideal Data and phase signal 6 waveform 


e 

OUTPUT 
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Binary ideal data is used throughout the simulations. For the ideal data, the logic 0 
and 1 have the same probability , it has a perfect symmetry(the duration of a digital one is 
equal to the time duration of a digital zero). Figure 3 shows the NRZ-L ideal binary data 
and waveforms of phase 6 (in radians) 



Figure 3 Phase Signals 6 Waveform 


Modulation 

The phase signal shown in Figure 3 can be modulated to an 8 PSK signal, s(t). To 
realize spectrum shaping, the phase signal, 6 , goes through a pulse Shaping filter before 
being modulated(type B). In this project, Bessel, Butterworth, and Square Root Raised 
Cosine filters are used. In addition, the change of the position of pulse shaping filters can 
produce a change of simulation results. If the filter is put behind the modulation[5], type A 
is simulated. This report will make a comparison to these two types of pulse shaping 
methods. Figure 4 shows the modulation scheme for unfiltered and filtered 8 PSK signals. 
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Phase 6 


modulator 
* A c e^ 


8PSK s(t) 

► 


(a) Unfiltered 8PSK With Constant Envelope 



(b) Filtered 8PSK Type A With Non-Constant Envelope 



(c) Filtered 8PSK Type B With Constant Envelope A c 

Figure 4. Modulation schemes for unfiltered and filtered phase signals, 
(b) is used in this project simulation. 


Pulse Shaping Filters 

Pulse Shaping Filters are used to narrow bandwidth and improve bandwidth 
utilization. In this report several types of filters such as 5th order Butterworth, 3rd Order 
Bessel and Square Root Raised Cosine (SRRC)(a=l) are utilized in the simulations for 
the ideal data. 

Pulse shaping can induce distortions of 8 PSK signals and increase the risk of Inter- 
Symbol Interference (ISI)[4]. The distortions due to the pulse shaping filters make the 
design of an optimal receiver difficult. 

Butterworth filters: Butterworth filters provide the maximally flat response at 
Q=0 and Q = co. Response varies monotonically during the stopband, decreasing smoothly 
from Q =0 and Q - oo . In the simulation, The 5th Order Butterworth filter is used. 

Bessel Filters: The Bessel type of filter has the slowest transition between 
passband and stopband and the passband phase response is better than the Butterworth. It 
has approximately linear phase in the passband so that it has a property of preventing 
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dispersion of the signal. The amplitude for Bessel filters is not as flat in the passband 
region as for the Butterworth filters. This project uses a 3rd Order Bessel Filter for the 
simulation. 

Square Root Raised Cosine Filter ( a =1) The Square Root Raised cosine filter 
obtains the filter coefficients according to the following formula[l]: 


[cos((1+a);r — )+ sin((1 -a)n —)] 
h(t)=4 a T — — T— 

K-J?[(4a U 2 - 1 ] 

In this equation, a is the roll off parameter and T is one symbol period(l /symbol rate). 
The Real Raised cosine block in SPW is selected for the project simulation. As mentioned 
in [5], if the Raised cosine filters are used the ISI will be eliminated in a linear channel. 
The SRRC filtering is noncausal and has infinite coefficients, but actually a finite number 
of coefficients are needed to determine the FIR tap length. This truncation can cause some 
frequency distortion. On the other hand, the distortion can be reduced by applying any of 
several windowing functions such as Bartlett, Blackman, or Hamming. However, filtering 
occurs before modulation, the channel is no longer linear as far as the filter is concerned. 

The Observation Of Unfiltered And Filtered 8 PSK Signals 


With the use of pulse shaping filters, the 8 PSK signals change their shape. The 
main character of the waveform change is a transient(rise and drop) in the filtered 6 (as 
shown in Figure 5). The transient line(rise or drop) makes a big difference between 
unfiltered(as shown in Figure 6) and the filtered 8PSK shown in Figure 7. 


unfiltered phase theta 




filtered phase theta 


10 0 0 
s a m p le s 



Figure 5 The Waveforms Of Unfiltered 6 And Filtered 6 
(Butterworth 5th filter is used, BT-1) 
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Figure 6 Real Part And Imaginary Part Of Unfiltered 8 PSK For Ideal Data 


real part of si ( t ) 




Figure 7 Real Part And Imaginary Part Of Filtered 8 PSK Type B (512 Samples/Symbol) 

From Figures 5 and Figure 7, it can be seen that the filtered phase 6 has a transient 
from 0 to 2 n . In sample range 0 to 500, the waveform of the filtered 8 PSK has a 
tremendous change in that transient range. It changes from a constant to a sinusoid wave. 
The reason for this change is the Intersymbol Interference(ISI) introduced by the filtering 
process. The filtering process elongates data symbols such that they overlap one another. 
When the pulse waveforms pass through the bandlimited system, they will be spread or 
dispersed. This spreading can cause ISI which increases error rates and affects power 
spectrum. 

Solid State Power Amplifiers SPA) 

A Solid State Power Amplifier (SSPA) is used in the project. The magnitude and 
phase plots of the SSPA were obtained from JPL. From the plot(Figure 8), the SSPA has 
two main ranges which include a linear and non-linear region. If the input power is less 
than -5 dB, the SSPA is worked in the linear region, and the phase offset is approximately 
zero. If the input power is larger than -5 dB, the SSPA is working in the nonlinear region. 
This project sets the SSPA to work at saturation level of 0 dB. 
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SSPA, ESA 10 Watts Amplifier: Magnitude Response Curve 




Figure 8 The Output Characteristic Of 10 Watts SSPA 


The main part of the SSPA is an analytical Traveling Wave Tube (TWT) block 
which can simulate the nonlinear channel. By looking up the SSPA table file provided by 
JPL, the TWT realizes a conversion from input amplitude to output amplitude and input 
amplitude to out phase. Figure 9 shows the structure of the SSPA. 
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Figure 9 The Structure of SSPA. 


Power Spectra Analysis 

Figure 10 shows the procedure for obtaining the power spectra for 8 PSK signal. 



Figure 10 Block Diagram of Spectrum Analysis 
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The random data source provides ideal data with a data rate of 256 Hz, and the 
sample frequency, Fs, is 131072 Hz. For the spectrum shaping, three types are used. They 
are 5th order Butterworth with BT=1, 2, 2.8, 3, 3rd order Bessel with BT=1, 1.2, 2, 3, and 
square-Root Raised Cosine with roll off factor a =1. The BT represents the product of the 
bandwidth with the symbol time. The constant envelope 8 PSK signal which comes from 
the modulator goes to the input of the Solid State Power Amplifiers SPA). 

The spectrum Analyzer can be put before and after the SSPA. Since the SSPA is 
worked at the saturation level and constant envelope 8PSK is used, the spectrum of the 
signal at the input of the SSPA is almost the same as that of the signal output from the 
SSPA. The Spectrum Analyzer is picked up from Interactive Simulation Library(ISL) of 
SPW. The power spectrum is obtained by taking the Fast Fourier Transform(FFT)of the 
input signals of the Spectrum Analyzer. The FFT length is 131072, and the window is 
Bartlett. The power spectrum is produced by taking 100 FFTs averaged together. The plot 
on the ISL window can be copied to the Signal calculator window by using Ctr-D. It is 
worthwhile to mention that the scaling on the vertical axis is relative. Suppose the vertical 
axis on the window is from 0 to -100 dB and the scale on the signal calculator window is -1 
to 1, one can obtain the relationship by using following formula 

S2=(S 1 - 1 )* 5 0 

Where S 1 is the value on the signal calculator window, and S2 is the actual value. 


Additive White Gaussian Noise Channel 


Additive White Gaussian Noise(AWGN) is added to the signals out of the SSPA to 
simulate channel noise. Figure 1 1 shows the power spectral density for the AWGN used in 
the simulations. Actual white noise has a flat spectrum from to °° . 

f PSD n (f) 


N 0 /2 



-f s /2 f s /2 

Figure 1 1 Power Spectral Density Of White Noise. 


9 

The AWGN is taken to have zero mean and with variance a n =N 0 /2. No/2 is the power 
spectral density of the white noise, and F s is the sampling frequency. 
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Receiver And Error Rate Estimation System 


There are three parts in the receiver and estimation system. These are the 
receiver(matched filter), synchronizer and error rate estimator. Figure 13 shows the block 
diagram for this system. 



Figure 12 Block Diagram For The Receiver Error Rate Estimator 


Receiver 

A matched filter is ideally used as a receiver. The matched filter is a linear filter 
which minimizes the noise and maximizes the signal. In this project, an integrate and dump 
block is used as the matched filter(matching to the NR Z data). Since the actual signal is 
fairly complex( due to the pulse shaping), a truly matched filter was not attempted at this 
time. 

Synchronizer 

The synchronizer performs a correlation between the reference signals and the 
received signals(delayed and distorted), and gives an estimate of the number of samples of 
delay and the phase between the two samples. The reference signal is then delayed the 
same number of samples so that the two signals can be compared to get the error rate. 
Figure 13 shows the block diagram for the synchronization. The DELAY & PHASE 
METER can detect the number of sample delays of the received signals. The COMPLEX 
VARIABLE DELAY lets the reference signals be delayed by the same number of samples 
as the received signal so they can be compared in the error rate estimator. 
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received 
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+ x Delay y 
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DELAY & PHASE 
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phase 

-► 
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+ delayed 

ready 
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Figure 13 Block Diagram of Synchronizer 


Error Rate Estimator 

The reference signal and received signal now having the same delay, are put into 
the error rate estimator to count the number of bit errors. The main block used to count bit 
errors is SIMPLE ERROR RATE ESTIMATOR(Figure 14)which compares the reference 
and the received signals to find bit errors. 


+ in SIMPLE rate 
ERROR 

+ ref RATE errs 
ESTIMATOR 


reset 

T 

Figure 14 Simple Error Rate Estimator Block Diagram 
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SIMULATIONS 


Figure 15 shows the system block diagram for the 8 PSK end to end system 
simulations. The simulations include two main parts: Power Spectrum and Bit Error Rates. 



Figure 15 Simulation System Block Diagram 


Simulations on Power Spectrum 

Figure 10 showed the location of where to obtain the power spectrum. The 
following parameters were used to obtain the power spectra. 

Data Source: 

a. Ideal data with p(0)( probability of zero )=0.5; 

b. Sample Rate (f s ) = 131072 Hz 

c. Symbol rate R s =256 (Bit rate=768; 512 samples/symbol); 

Base band Filters 

a. None 

b. Butterworth 5th Order (BT=1,2,2.8,3) 
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c. Bessel 3rd Order (BT=1,1.2,2,3) 

d. SRRC(« =1, Number of tap length- 16* 5 12=8 192, Bartlett window, 

Square-root raised cosine=‘yes’; Sample input=‘yes’) 

Power Amplifier: 

European Space Agency(ESA) 10-watt SSPA 

SPW Spectrum Analyzer 

a. FFT length=131072 

b. Window type: Bartlett 

c. Average: 100 FFTs 

d. Sampling frequency =131072 

Spectrum Plots ranges 

±10 Rt,; ± 75 Rb; ± 250 Rb 


The plots of the power spectra are given in Figures 16-34. 
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SSPA Output, +/-10Rb (5th Order Butterworth BT=1, Ideal Data) 



(a) 5th Order Butterworth Filter (BTM), ± 10 Rb 


SSPA Output, +/-75Rb (5th Order Butterworth BT=1 , Ideal Data) 



(b) 5th Order Butterworth Filter (BTM) , ± 75Rb 


SSPA Output, +/-250Rb (5th Order Butterworth BT=1, Ideal Data) 



(c) 5th Order Butterworth Filter (BT=1), ± 250Rb 
Figure 17: 5 th Order Butterworth Filter (BT=1) Power Spectra, With Central Spike 
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SSPA Output, +/-10Rb (5th Order Butterworth BT=1, Ideal Data) 
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-40 
co -60 
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-75 -60 -40 -20 0 20 40 60 75 
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( b) 5th Order Butterworth Filter (BT=1) , ± 75Rb 
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g -60 
-80 
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Frequency (Rb),Rb=256*3, Without Max Spikes 


SSPA Output, +/-250Rb (5th Order Butterworth BT=1 , Ideal Data) 
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i i i 1 
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(a) 5th Order Butterworth Filter (BT=1), ± 10 R b 

SSPA Output, +/-75Rb (5th Order Butterworth BT=1, Ideal Data) 




(c) 5th Order Butterworth Filter ,BT=1, ± 250R b 
Figure 18: 5 th Order Butterworth Filter Power Spectra, Without Central Spike 
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SSPA Output, +/-10Rb (5th Order Butterworth BT=2, Ideal Data) 



(a) 5th Order Butterworth Filter (BT-2) , ± 10Rb 


SSPA Output, +/-75Rb (5th Order Butterworth BT=2, Ideal Data) 



SSPA Output, +/-250Rb (5th Order Butterworth BT=2, Ideal Data) 



(c) 5th Order Butterworth Filter (BT=2), ± 250Rb 
Figure 19: 5 th Order Butterworth Filter (BT=2) Power Spectra, With Central Spike 
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SSPA Output, +/-10Rb (5th Order Butterworth BT=2, Ideal Data) 



(a) 5th Order Butterworth Filter (BT-2) , ± 10Rb 


0 

-20 

-40 

§ -60 
-80 
-100 


-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3, Without Max Spikes 

(b) 5th Order Butterworth Filter (BT=2) , ± 75Rb 


SSPA Output, +/-75Rb (5th Order Butterworth BT=2, Ideal Data) 



SSPA Output, +/-250Rb (5th Order Butterworth BT=2, Ideal Data) 



(c) 5th Order Butterworth Filter (BT=2), ± 250Rb 


Figure 20: 5 th Order Butterworth Filter (BT^), Without Central Spike 
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(a) 5th Order Butterworth Filter (BT=2.8), ± 10Rb 

SSPA Output, +/-75Rb (5th Order Butterworth BT=2.8, Ideal Data) 



Frequency (Rb),Rb=256*3,With All Spikes 

(b)5th Order Butterworth Filter (BT=2.8), ± 75Rb 
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Frequency (Rb),Rb=256*3,With All Spikes 


(c) 5th Order Butterworth Filter (BT-2.8), ± 250Rb 
Figure 21 : 5 th Order Butterworth Filter (BT=2.8) Power Spectra, With Central Spike 


27 






SSPA Output, +/-10Rb (5th Order Butterworth BT=2.8, Ideal Data) 



(a)5th Order Butterworth Filter (BT=2.8), ± 10Rb 


SSPA Output, +/-75Rb (5th Order Butterworth BT=2.8, Ideal Data) 



Frequency (Rb),Rb=256*3, Without Max Spikes 

(b)5th Order Butterworth Filter (BT=2.8), ± 75Rb 


SSPA Output, +/-250Rb (5th Order Butterworth BT=2.8, Ideal Data) 
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-80 

-100 


-250 -200 -150 -100 -50 0 50 100 150 200 250 

Frequency (Rb),Rb=256*3, Without Max Spikes 

(c) 5th Order Butterworth Filter, BT=2.8, ± 250Rb 
Figure 22: 5 th Order Butterworth Filter, (BT-2.8) Without Central Spike 
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Frequency (Rb),Rb=256*3,With All Spikes 

(a) 5th Order Butterworth Filter (BT-3), ± 10Rb 



SSPA Output, +/-75Rb (5th Order Butterworth BT=3, Ideal Data) 



(b) 5th Order Butterworth Filter (BT=3), ± 75Rb 


SSPA Output, +/-250Rb (5th Order Butterworth BT=3, ideal Data) 



(c) 5th Order Butterworth Filter (BT=3), ± 250Rb 
Figure 23: 5 th Order Butterworth Filter (BT=3) Power Spectra, With Central Spike 
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SSPA Output, +/-10Rb (5th Order Butterworth BT=3, Ideal Data) 



(a) 5th Order Butterworth Filter (BT-3), ± 10 Rb 
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(b) 5th Order Butterworth Filter (BT-3), ± 75Rb 
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SSPA Output, +/-250Rb (5th Order Butterworth BT=3, Ideal Data) 
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(c) 5th Order Butterworth Filter (BT-3), ± 250Rb 
Figure 24: 5 th Order Butterworth Filter (BT=3) Without Central Spike 
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SSPA Output , +/-10Rb (3rd Order Bessel BT=1, Ideal Data) 



Frequency (Rb),Rb=256*3,With All Spikes 

(a) 3rd Order Bessel Filter (BT=1) , ± 10Rb 


SSPA Output , +/-75b (3rd Order Bessel BT=1 ,ldeal Data) 


oo 

■d 


-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3,With All Spikes 



(b) 3rd Order Bessel Filter (BT-1), ± 75Rb 


SSPA Output, +/-250Rb (3rd Order Bessel BT=1, Ideal Data) 



-250 -200 -150 -100 -50 0 50 100 150 200 250 

Frequency (Rb),Rb=256*3,With All Spikes 


(c) 3rd Order Bessel Filter (BT— 1) , ± 250Rb 
Figure 25: 3 rd Order Bessel (BT=1) Power Spectra, With Central Spike 
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0 

-20 

-40 


SSPA Output, +/-10Rb, 3rd Order Bessel BT=1 , Ideal Data) 


-10 -8 -6 -4 -2 0 2 4 6 8 10 

Frequency (Rb),Rb=256*3, Without Max Spikes 

(a) 3rd Order Bessel Filter (BT=1) 9 ±10Rb 





(b) 3rd Order Bessel Filter (BTM), ± 75 Rb 


SSPA Output, +/-250Rb, 3rd Order Bessel BT=1, Ideal Data) 



-250 -200 -150 -100 -50 0 50 100 150 200 250 

Frequency (Rb),Rb=256*3, Without Max Spikes 


(c) 3rd Order Bessel Filter (BT=1), ±250Rb 
Figure 26: 3 rd Order Bessel Filter (BT=1), Without Central Spike 
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(b) 3rd Order Bessel Filter (BT-1.2) , ± 75Rb 


SSPA Output, +/-250Rb (3rd Order Bessel BT=1.2,ldeal Data) 



(c) 3rd Order Bessel Filter (BT-1.2), ± 250Rb 
Figure 27: 3 rd Order Bessel Filter (BT=1.2) With Central Spike 
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-80 

-100 


-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3, Without Max Spikes 

(b) 3rd Order Bessel Filter (BT=1.2), ± 75Rb 


SSPA Output, +/-10Rb, 3rd Order Bessel BT=1.2,ldeal Data) 



(c) 3rd Order Bessel Filter (BT=1.2), ± 250Rb 
Figure 28: 3 rd Order Bessel Filter (BT=1.2), Without Central Spike 
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SSPA Output , +/-10b (3rd Order Bessel BT=2, Ideal Data) 



(a) 3rd Order Bessel Filter (BT-2), ± 10Rb 



j i i -na—^ i i 1 1 

-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3,With All Spikes 


(b) 3rd Order Bessel Filter (BT-2) , ± 75Rb 


SSPA Output, +/-250Rb (3rd Order Bessel BT=2, Ideal Data) 



(c) 3rd Order Bessel Filter (BT-2), ± 250Rb 
Figure 29: 3 rd Order Bessel Filter (BT=2), With Central Spike 


35 






-80 
-100 

-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3, Without Max Spikes 

(b) 3rd Order Bessel Filter (BT=2), ± 75Rb 



-250 -200 -150 -100 -50 0 50 100 150 200 250 

Frequency (Rb),Rb=256*3, Without Max Spikes 


(c) 3rd Order Bessel Filter (BT=2), ± 250Rb 
Figure 30: 3 rd Order Bessel Filter (BT=2) Without Central Spike 
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SSPA Output, +/-10Rb (3rd Order Bessel BT=3, Ideal Data) 



(a) 3rd Order Bessel Filter (BT=3), ± 10Rb 


SSPA Output, +/-75Rb (3rd Order Bessel BT=3, Ideal Data) 



-75 -60 -40 -20 0 20 40 60 75 

Frequency (Rb),Rb=256*3,With All Spikes 


(b) 3rd Order Bessel Filter (BT-3), ± 75Rb 



-250 -200 -150 -100 -50 0 50 100 150 200 250 

Frequency (Rb),Rb=256*3,With All Spikes 


(c) 3rd Order Bessel Filter (BT=3), ± 250Rb 
Figure 31: 3 rd Order Bessel Filter (BT=3), With Central Spike 
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SSPA Output, +/-10Rb (3rd Order Bessel BT=3, Ideal Data) 



Frequency (Rb),Rb=256*3, Without Max Spikes 

(a) 3rd Order Bessel Filter (BT=3), ± 10 Rb 


SSPA Output, +/-75Rb (3rd Order Bessel BT=3, Ideal Data) 



(b) 3rd Order Bessel Filter (BT=3),. ± 75Rb 


SSPA Output, +/-250Rb (3rd Order Bessel BT=3, Ideal Data) 



(c) 3rd Order Bessel Filter (BT=3), ± 250Rb 
Figure 32: 3 rd Order Bessel Filter (BT=3) Without Central Spike 






SSPA Output, +/-10Rb (Raised Cosine, a=1, Ideal Data) 



SSPA Output, +/-75Rb (Raised Cosine, a=1, Ideal Data) 



(b) Square-Root Raised Cosine ( a -1) , ± 75Rb 


SSPA Output, +/-250Rb (Raised Cosine, a=1, Ideal Data) 



(c) Square-Root Raised Cosine ( a =1), ± 250Rb 
Figure 33: Square-Root Raised Cosine (a =1) With Central Spike 
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SSPA Output, +/-10Rb (Raised Cosine, a=1, Ideal Data) 



(a) Square-root Raised Cosine ( a =1) , ± 10 Rb 


SSPA Output, +/-75Rb (Raised Cosine, a=1, Ideal Data) 



(b) Square-root Raised Cosine (a =1), ± 75Rb 


SSPA Output, +/-250Rb (Raised Cosine, a=1 , Ideal Data) 



(c) Square-Root Raised Cosine ( a =1) , ± 250Rb 
Figure 34: Square-Root Raised Cosine (a =1), Without Central Spike 
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Frequency Band Utilization Ratio ( p ) 


The frequency band utilization ratio can be defined as: 

p- Number of spacecraft with Filtering Accommodated in Frequency Band 
Number of spacecraft without Filtering Accommodated in Frequency Band 

From the simulations, the utilization ratios (type B) were calculated and listed in 
Table 2 along with the utilization ratios for Type A[5]. From the utilization ratio table, for 
the same BT and baseband filters, modulator type B has a higher frequency band 
utilization ratio than type A. 


Utilization Ratio (Type A) 


Filter Type 

-50 dB point 

utilization 

None 

35Rb 

1 

Butterworth(5th order)BT=T 

2.5Rb 

14 

Beset (3rd order) BT=1 

2.95Rb 

11.86 

SRRC rolloff=l 

2.5Rb 

14 


Utilization Ratio (Type B) 


Filter Type 

-50 dB point 

utilization 


with spikes 
(Rb) 

no spikes 

(Rb) 

with spikes 

no 

spikes 

None 

35 

1 

1 

Butterworth(5th order) BT=1 

1.22 

1.28 

28.7 

27 

Butterworth(5th order) BT=2 

2.38 

2.43 

14.7 

14.4 

Butterworth (5th order) BT=2.8 

3.2 

3.23 

11 

10.8 

Butterworth (5th order) BT=3 

3.45 

3.45 1 

10.1 

10.1 

Beset (3rd order) BT=1 

1.625 

1.75 

21.5 

20 

Beset (3rd order) BT=1.2 

1.88 

2.04 

18.6 

17.1 

Beset (3rd order) BT=2 

3.1 

3.22 

11.3 

10.8 

Beset (3rd order) BT=3 

4.5 

4.53 

7.8 

7.7 

SRRC rollofW 

1.11 

1.23 

31.5 

28.5 


Table 2 Utilization Ratio of Type A and B 8 PSK 
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Bit Errors Rates Simulations 


The bit errors are produced due to the ISI and AWGN. This report provides the 
measure of bit errors according to bit energy to noise ratio(Eb/N 0 ).Figure 15 shows the 
system block diagram for measuring bit error rates. For different filter types and their 
positions(type A or B), the BERs are put together for comparison. The BER plots are 
given in Figure 35-40. 

The following parameters were used to obtain the bit error rates: 

Data source: 

a. Ideal data with p(0)( probability of zero )=0.5; 

b. Sample Rate (f s ) = 16 Hz 

c. Symbol rate R s =l (Bitrate=3; 16 samples/symbol); 

Base band Filters 

a. Butterworth 5th Order (BT=1,2,2.8,3) 

b. Bessel 3rd Order (BT=T,1.2,2,3) 

c. SRRC( a =1, Number of tap length^ 16* 16=256, rectangular 

window. Square-root raised cosine^ y es ^ Sample input=‘yes’) 


Power Amplifier: 

European Space Agency(ESA) 10-watt SSPA , Back off =0dB as obtained 
from JPL 



Figure 35. BER Plot For 5th Order Butterworth Filters(Type B) 
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Bit Error Probability 
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Bit Error Probability Bit Error probability 


BER ( 5th Order Butterworth Filter; SSPA in Saturation Level) 



Figure 38. BER Plot For 5th Order Butterworth Filters(T.ype A And B) 


BER ( 3rd Order Bessel Fi!ter;SSPA in Saturation Level) 



Eb/No(dB) 

Figure 39 BER For 3rd Order Bessel Filters (Type A And B) 
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BER (SRRC Filter;SSPA in Saturation Level) 



Figure 40 BER Plots For SRRC Filters (Type A And B) 


ISI Losses Due To Filter 


The ISI losses due to the pulse shaping filter can be determined by the BER 
simulations. The following tables gives the results of baseband filter ISI losses at 10° BER 
compared to theoretical 8 PSK(or BPSK) and the filtered 8 PSK. 


Filter Type 

Loss 

BT=1 

(dB) 

Loss 

BT=1.2 

(dB) 

Loss 

BT=2 

(dB) 

Loss 

BT=2.8 

(dB) 

Loss 

BT=3 

(dB) 

Butterworth 
5th Order 

22.5-10= 

12.5 


12.32-10= 

2.32 

11.48-10= 

1.48 

11.48-10= 

1.48 

Bessel 
3rd Order 

1 

l 

t 

I 

l 

1 

l 

1 

l 

l 


12.85-10= 

2.85 


11.20-10= 

1.20 

SRRQ>=1) 

14.22-10= 

4.22 










Table 3 Baseband Filters ISI Loss Measurements at 10' 3 BER(Compared to ideal 8PSK) 
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Filter Type 

Loss 

BT=1 

(dB) 

Loss 

BT=1.2 

(dB) 

Loss 

BT=2 

(dB) 

Loss 

BT=2.8 

(dB) 

Loss 

BT=3 

(dB) 

Butterworth 
5th Order 

22.5-6.8= 

15.7 


12.32-6.8= 

5.52 

11.48-6.8= 

4.68 

11.48-6.8= 

4.68 

Bessel 
3rd Order 

K 

1 

1 

1 

I 

I 

1 

1 

1 

1 

20-6.8= 

13.2 

12.85-6.8= 

6.05 


11.20-6.8= 

4.88 

SRRQ>=1) 

14.22-6.8= 

7.42 










Table 4 Baseband Filters ISI Loss Measurements at 10‘ 3 BER(Compared to ideal BPSK) 
(dashes indicate that the simulation was not performed) 

As stipulated in [9], the optimum filter should have the smallest BT value and provide 
losses<0.4 dB. Thus a higher order of BT should be simulated to find an optimum BT and 
meet the requirement of ISI loss <0.4. The Figure 41 shows the relationship of ISI loss and 
utilization ratio for different filters. 



Utilization Ratios 

Figure 41. Plots For ISI Loss As Utilization Ratios (Type B) 

( for 5th Order Butterworth Filter ; for 3rd Order Bessel Filter; * for SRRC filter) 
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CONCLUSIONS AND SUGGESTIONS 

The simulations conducted in the Space Communication and Telemetry Laboratory 
at NMSU on the pulse shaped 8 PSK modulations were concluded as follows. 

Different types and positions of baseband filters have different effects for 
narrowing bandwidth or frequency band utilization ratio. In general, type B had better 
results for narrowing the bandwidth than type A did for the same BT value. 

Different types and positions of baseband filters also have different effects on 
producing bit error rates. In general, type A had lower bit error rates than type B. In 
general, that means that the narrower the bandwidth ,the larger the bit error rate. For the 
3rd Order Bessel filter, the BER produced by type A(BT=2) is almost the same as that 
produced by type B (BT=3). For the 5th Order Butterworth filters, the BER produced by 
type B(BT=2.8 or 3) is between those produced by type A (BT=1) and type A(BT=2). 

A matched filter is crucial for receiving the signal to lower the BER. Integrate and 
dump circuits can exactly match pulse signals, but for pulse shaped signals it could not 
match perfectly. A new matched filter which can match pulse shaped signals is needed so 
that the BER can be reduced further. More filter types and larger BT(BT>4) should be 
simulated to find an optimal filter and BT. The Spikes that are prevalent in the power 
spectra should be investigated and ways to reduce or eliminate them should be found. 

Work has been performed with 8 PSK (with no pulse-shaping) passed through a non-linear 
satellite channel consisting of a prefilter, Traveling Wave Tube Amplifier, (TWTA), and a 
post filter [14]. The pre and post filters were 6 th order Butterworth filters with a BT=0.75. 
Equalizer methods were investigated to counter the nonlinear ISI created by the system. 
These equalizer techniques included a fixed Volterra Equalizer adapted with an Least 
Mean Square (LMS) algorithm, an Adaptive Fractionally-Spaced Equalizer (FSE) with the 
Volterra Equalizer, a Single Filter (SF) Maximum-Liklihood Sequence Detector (MLSD), 
and a matched filter bank MLSD (MFB-MLSD). These same techniques could be used to 
counter the non-linear ISI created by the pulse shaping and therefore decrease the BER. A 
block diagram of the system and some of the results from [14] are given in Figure 42, 
however for a detailed description see [14]. 
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w(t) 



(a) Communication System Model Used for Equalization 



(b) Performance of MFB-MLSD receiver compared to FSE-Volterra receiver 



(c) Performance of SF-MLSD and MFB-MLSD receivers 


Figure 42: Some Results of Simulations Using Equalizers 
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